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Abstract. The description of many dynamical problems like the particle motion in higher dimensional 
spherically and axially symmetric space-times is reduced to the inversion of a holomorphic hyperellip- 
tic integral. The result of the inversion is defined only locally, and is done using the algebro-geometric 
techniques of the standard Jacobi inversion problem and the foregoing restriction to the 6-divisor. For a 
representation of the hyperelliptic functions the Klein- WeierstraB multivariable sigma function is intro- 
duced. It is shown that all parameters needed for the calculations like period matrices and Abclian images 
of branch points can be expressed in terms of the periods of holomorphic differentials and theta-constants. 
The cases of genus two and three are considered in detail. The method is exemplified by particle motion 
associated with a genus three hyperelliptic curve. 
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1. Introduction 

In a wide range of dynamical problems of classical systems one faces the problem of the inversion of 
integrals of the type |G02] 
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where V n is a polynomial of order n. Such so-called hyperelliptic integrals of genus g, n = 2g + 1, are 
found, for example, for the particle motion in higher dimensional axially symmetric space-times. The 
solution of the inversion is the function x = x(t). If the genus is g = 1 then the result of the inversion 
is an elliptic function, which is a doubly periodic function of one complex variable. For higher genera 
g > 1 such an inversion becomes impossible because, as already recognized by Jacobi (see e.g. [Mark92] ). 
2g-periodic functions of one variable do not exist. However, Jacobi was able to resolve this contradiction 
by formulating his celebrated Jacobi inversion problem that involves g hyperelliptic integrals. The problem 
was solved in terms of so-called hyperelliptic functions which are indeed 2g-periodic functions for g > 1 
which depend on g variables while the periods (also called moduli) are g x ^-matrices. The domain of 
these hyperelliptic functions - the Jacobi variety - is thus the (/-dimensional complex space C 9 factorized 
by the period lattice. 

The Jacobi inversion problem stimulated the development of algebraic geometry and in particular led 
to the discovery of solutions of many classical mechanical systems like Neumann's geodesic on an ellipsoid, 
the spinning top of Kowalewskaja, Kirchhoff's motion of a rigid body in a fluid, and others that were 
integrated using Jacobi's procedure. This special type of integrability that might be called algebro- 
geometric integrability has been receiving much attention owing to the discovery of the vast class of partial 
differential equations of Korteveg-de Vries type that admits this type of integrability. 

In this paper we consider the problem of the inversion of one hyperelliptic integral on the basis of 
the well developed algebro-geometric technique for the standard Jacobi inversion problem. The results of 
such inversions can be obtained by a restriction of hyperelliptic functions to special subsets of the Jacobi 
variety - the 9 -divisor, that is given as a solution of the equation including the Riemann 9 -function. Such 
restrictions of hyperelliptic functions can be defined only locally; alternatively they can be realized on an 
infinitely sheeted Riemann surface [FG07] . These functions inherit a number o f prope rties of standard 



elliptic functions, like the addition formulae of the Frobenius-Stickelberger type [Oni02 



Our development starts with the standard Jacobi inversion problem that involves g hyperelliptic inte- 
grals with variable bounds called divisor and describes a dynamic system with g degrees of freedom. Then 
we are fixing 1 < m < g points of the divisor making it special or a divisor with deficiency. In the context 
of this paper that means that we are considering the case when the genus of the underlying algebraic curve 
exceeds the number of degrees of freedom of the system. Although different values of m appear in various 
problems we are concentrating here on the case of maximal deficiency m = g — 1, i.e., on the inversion of 
one hyperelliptic integral. 

The Jacobi inversion problem for a divisor with deficiency has a long history that includes Baker's 
consideration Bak907 , Grant's [Gra90 and Jorgenson's [Jor92 treatment of the genus two case, Onishi's 
consideration 



Oni98 



of the genus three case, description of certain dynamic systems with separable vari- 
ables [EEKL93 , the treatment of the weak Kowalevski-Painleve property |AF00j , the integration of Somos 
sequences MatQ3., BEH05 , }Hon07j . reductions of Benney hierarchies BG04 , BG06] and others. Here 
we will consider the problem of inversion in a systematic way within the Klein- Weierstrafi realization of 
the theory of Abelian functions that is documented in the book of Baker [Bak 897] (see also the review 
[BEL97] and the more recent developments in Buchstaber and Leykin }BL05j . also Nakayashiki |Nak08a] . 
Matsutani and Previato jMPlOj ). The inversion formulae discussed here result in the restriction of the 
solution of the standard Jacobi inversion problem written in terms of the Kleinian p-functions to the corre- 
sponding stratum of the 0-divisor. In this context there appears the problem of a suitable parametrization 
of the 0-divisor in the case of higher genera. We solve this problem on the basis of Newton's method for 
the approximation of multivariable functions. We note that in this paper we are discussing the inversion 
of holomorphic integrals only; similar considerations can be undertaken for meromorphic integrals and 
integrals of the third kind, see e.g. |EPR03j . 

In order to elaborate an effective calculation procedure for the inversion of one holomorphic hyperelliptic 
integral we are solving a problem of general interest, that is, the calculation of the period matrix of 
meromorphic differentials that is sometimes called second period matrix. Existing Maple codes contain 
the evaluation of the Riemann period matrix by the given curve, i.e., period matrices of holomorphic 
differentials, only. However, the calculation of the periods of meromorphic differentials is unavoidable in 
certain problems and in particular in our inversion problem. Therefore we show that it is possible to express 
this second period matrix in terms of the first period matrix and theta-constants. In this way we reduce 
the number of complete integrals, which are necessary to calculate, to the Riemann period matrix given by 
Maple codes. Another result yields the characteristics of the Abelian images of branch points by the given 
holomorphic period matrices and, thus, that the homology basis can be reconstructed by these data. That 
permits to find the vector of Riemann constants and to carry out the whole calculation without referring 
to the homology basis. While ^-functional calculations are usually considered as technically complicated, 
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here we describe an easily algorithmized scheme that turn such calculations into routine procedures at 
least in the hyperelliptic case. 

Our paper is organized as follows. In Section [2] we shortly describe the problem of particle motion in 
higher dimensional spherically and axially symmetric space-times. This problem will serve as a laboratory 
for the approbation of the methods developed. In Section [3] we recall the known facts from the theory 
of hyperelliptic functions and develop a realization of these functions in terms of Klein- Weierstrafi multi- 
variable er-function. Special attention is focused on the effective calculation of the moduli of the system. 
We show in particular how to express periods of meromorphic differentials in terms of the theta-constants 
and periods of holomorphic differentials. In Section 0] we are considering the stratification of the ^-divisor 
and show how to single out the stratum that is the image of the curve inside the Jacobian in terms of con- 
ditions on the cr-function. This stratum serves as the domain for the quasi-elliptic function. In Sections [5] 
and [5] we are considering the application of the method developed to the cases of genus two and genus 
three hyperelliptic curves. Finally in the last Section[7]we are coming back to the initial physical problem 
and demonstrate how to compute the trajectories of test particles in higher dimensions by the method 
of restriction to the #-divisor. We explicitely calculate orbits in a 9-dimensional Reissner-Nordstrom-de 
Sitter space-time, which is characterized by its mass, electric charge and the cosmological constant. The 
underlying polynomial is of degree 7 which corresponds to a genus 3 hyperelliptic curve. This is the gen- 
eralization of the examples considered in HKKL08 , where orbits with underlying hyperelliptic curves of 
genus 2 where calculated. 

We believe that our method has much wider applications than the special physical problem considered 
here and that it can be used in other problems that needs the inversion of a hyperelliptic integral. The 
same approach works for the inversion of meromorphic integrals, that we will consider elsewhere. Some of 
our results can be used for the Jacobi inversion problem on the strata with smaller deficiency. 

2. Particle motion in General Relativity 

Ordinary differential equations of the form 

d:r 



(2-1) — = f{x,y/P&)), 

where / is a rational function of x and V^(x) is a polynomial of order four are solved by elliptic integrals 
introduced by Jacobi and WeierstrafQ. The corresponding equations of that form with a fourth order 
polynomial Pi(x) can be reduced to one with a third order polynomial. As an example we mention the 
motion of a point particle or light ray given by the geodesic equation 

, _ d x^ r ii dx^ dx „ 

with the Christoffel symbol 

(2-3) { ^ } := -g^ {d p g al , + d a g pv - d v g pa ) , 

where g^ is the space-time metric, ds is the proper time defined by ds 2 = g^dx^dx 11 . For a light-like 
particle (photon) the parameter s is replaced by some affine parameter. 
For a Schwarzschild metric 

(2.4) ds 2 = g tt dt 2 - g rr dr 2 - r 2 (dd 2 + sin 2 Mip 2 ) 

with 

, . 1 2M 
2.5 g tt = — = 1 , 

g rr r 

where M is the mass of the gravitating body (we choose units so that the Newtonian gravitational constant 
as well as the velocity of light are unity, G = c = 1), this geodesic equation with a substitution x = f(r) 
yields an equation describing the dependence of the radial coordinate r on the azimuthal angle ip 

(2.6) G|) 2= 4x3 - 92X ~ 93 > 

where the Weierstrafi invariants g2 and g% depend on M and the energy and angular momentum of 
the particle. The complete set of solutions in terms of the Weierstrafi p-function have been given and 
extensively discussed by Hagihara |Hag31[ . The corresponding periods of the p-functions are directly 



^Here and below we will closely follow the standard notations of the theory of elliptic functions fixed in BE55 . 



4 



V.Z. ENOLSKI, E. HACKMANN, V. KAGRAMANOVA, J. KUNZ, AND C. LAMMERZAHL 



related to observable effects like the perihelion shift and the deflection (scattering) angle of massive bodies 
and of light. 

Geodesic equations of neutral test particles or photons in Taub-NUT space-times [KKHL10 and of 
charged test particles in Reissner-Nordstrom space-times [GKlOl are also solved in terms of elliptic Weier- 
strafi functions. In the case of a Schwarzschild-(anti-)de Sitter or Kerr-(anti-)dc Sitter metric we encounter 
similar equations but with a polynomial of fifth and sixth order (where the sixth order polynomial can be 
reduced to a fifth order one) 

(2.7) (x^ = V 5 (x) ,i = 0,l, or, f^j = (x - c) 2 V 5 (x) , 

where c is a constant. The corresponding equations have been solved explicitely in [HL08 ( HL08a, HKKL09, 
IHKKL09aj . 

These examples can be generalized further to include polynomials of even higher orders as outlined in 
the following. 

2.1. Geodesic equations in higher dimensional spherically symmetric space times. The metric 
of a spherically symmetric Reissner-Nordstr6m-(anti-)de Sitter space-times in d dimensions is given by 

(2.8) ds 2 = g tt dt 2 - g rr dr 2 - r 2 dQ 2 d _ 2 , 
with 

1 fr s \d-3 2Ar 2 /o\2(<i-3) 

< 2 - 9 > ^ = ^ = 1 -(f) - (rf -f) (rf -2) + (;) - 

where A is the cosmological constant, q the charge of the gravitating mass M, r$ = 2M, and dVl 2 d _ 2 is the 
surface element of the d — 2-dimensional unit sphere. The geodesic equation then leads to 

„ 10) 

\d(pj L z g rr g tt \ \ 



r 4 / 2 / /r s \ d -3 2Ar 2 /q\ 2 ( rf -3)\ f „ L 2 



(2 - U) ~£ 2 V V ^rJ (d-l)(d-2) ' \r) J V ' r 2 

where E and L are two conserved quantities: the dimensionless energy E, and the angular momentum L 
with the dimension of length (both are normalized to the mass of a test particle) 

L ^t; 

where A is an affine parameter along the geodesic. 8=1 for massive test particles and (5 = for massless 
particles. A substitution x = f(r) gives equations of the form 

; dx x 2 



(2.13) \fd^) =Vni - x) 

for some < i < g where V n denotes a polynomial of order n and g = f 2 ^] is the genus of a curve 
w 2 = V n {x). In some cases through appropriate substitutions the order of the polynomial can be reduced 
HKKL08 . However, in general we have n > 7, as in the example in Section [7] 

2.2. The effective one— body problem. Another example of spherically symmetric problems is the rela- 
tivistic effective one-body problem in four dimensions. While in Newtonian gravity the two-body problem 
can be exactly reduced to an one-body problem this is not possible in Einstein's General Relativity. The 
rclativistic two-body problem can be reduced to a one-body problem only in terms of a series expansion. 
In this framework the relative coordinate between two bodies with masses Mi and M 2 formally fulfills a 
geodesic equation in a space-time with the effective metric 

(2.14) ds 2 = -g tt (r, v)dt 2 + g rr (r, v)dr 2 + r 2 (dd 2 + sin 2 Mip 2 ) , 
where u = 2(M 1 + M 2 )/r, v = MiM 2 /(Mi + M 2 ) 2 and 

g tt (r, v) = 1 - 2u + 2vv? + va^ + <&(u 5 ) 
^ 2 ' 15 ^ (gtt(r, v)g rr {r 1 v))' 1 = 1 + Qvu 2 + 2(26 - 3is)isu 3 + . 

The corresponding effective one-body equation of motion is then given by [BD991 [DJS00 

dr\ 2 r 4 1 (. L 2 



(2.16) — =— \E 2 ~g tt [l + — 

dtp J L z g rr gtt V V r 
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where E and L are again the conserved energy and angular momentum. This is a series expansion which 
can be expanded to arbitrary order. Although this is only a series expansion, analytic methods are helpful 
for the purpose to have a complete discussion of the possible types of orbits of a binary system. 

2.3. Geodesic equations in higher dimensional axially symmetric space times. As an exam- 
ple of d-dimensional axially symmetric space-times we consider the simplest one, namely the Myers- 
Perry [MP86] space-times with only one rotation parameter a given by |KKZ10( |Emp08| 

1 / 2M A , 2 4aMsin 2 6> , , sin 2 6 ( , 9 2x 2 2a 2 M . 2 \ , , 



ds z =- 



( 2.17) + ^dr 2 + p 2 dQ 2 + r 2 cos 2 6dn 2 n , 



where p 2 = r 2 + a 2 cos 2 8, A = (r 2 + a 2 ) — Jj£f ± , and n — d — 4. Here, M is the mass of the black hole and 

the surface element of the unit n-sphere is df2 2 = Yn=i 111=1 sm2 Tpkdip 2 . 

Owing to the conservation laws related to the symmetries of the underlying space-time there are a 
conserved energy E, an angular momentum and n further constants ty 2 , i — 1, . . . , n. As a consequence, 
it is possible to separate the Hamilton- Jacobi equation [VSP05] . The resulting equations of motion are 
then given by 

(2.18) = Vfl 

(2.1 9 ) = Ve 

<"°» fx = i^+^-^Tm+^r^-^^TsTx 



dt r 2 + a 2 r , o 1 dr r , „ , 1 dfl 

(2.21) - = —[ir 2 + a 2 ) E-aL v ]^- + a[L^a^ 2 eE] 7 =- 

, n „„s dibi \J~A~i 1 

dA ffi=i sm 2 V'fc r cos 61 



with A, = * 2 - and 

(2.23) JJ(r) = [(r 2 + a 2 )£ - ai v ] 2 - A (k + 5r 2 + ^ 2 ^j , 

(2.24) 9(0) = K - Sa 2 cos 2 9 - -^f- %- \L V - a sin 2 BE] 2 , 

cos^ sin 6 

where K is a further constant, called the Carter constant, which emerges from the separation process. It 
is obvious that V(r) in R(r) — Vd r t,-^ is a polynomial whose order increases with the dimension of the 
space-time. (In some cases the order of V{r) can be reduced by a substitution.) 

In the following sections we will explain the theory and the details of how to analytically solve the 
geodesic equations in the above cases and for similar physical problems. 



3. HYPERELLIPTIC FUNCTIONS 

The solutions of the differential equations given in the foregoing section can be considered as points on 
a hyperelliptic curve X g of genus g given by the equation 

2g+l 2g+l 

(3.1) w 2 = v 2g+ i (z) = x ^ = 4 n ( z - efc )' 

i=0 fc=l 

and realized as a two sheeted covering over the Riemann sphere branched in the points (efc,0), k 6 Q = 
{1, . . . , 2(7+1}, with ej ^ for j ^ k, and at infinity, e 2g +2 = oo. Notice that we do not require the eu to be 
real. However, when they are real, we find it convenient to order them according to e± < e 2 < ■ ■ ■ < e 2g +i, 
i.e., in the opposite way as compared to the Weierstrafi ordering, see Fig. [TJ Denote P = (z,w) a 
coordinate of the curve. The factor 4 in (|3.1[) is introduced to preserve resemblance with the Weierstrafi 
cubic for g = 1, 

(3.2) w 2 = 4z 3 - g 2 z - g 3 = 4(« - ei)(« - e 2 )(z - e 3 ). 
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Figure 1 . A homology basis on a Riemann surface of the hyperelliptic curve of genus g 
with real branch points e%, . . . , &2g+2 — oo (upper sheet). The cuts are drawn from e%i— 1 
to e2i for i = 1, . . . , g + 1. The b-cycles are completed on the lower sheet (the picture on 
the lower sheet is just nipped horizontally). 

As shown in Fig. Q] we equip the hyperelliptic curve X g with a canonical homology basis 
(3.3) (oi, ...,a g ;bx,...,b g ), cu o b; = -bi o a ( = <% , en o aj = b, o b, = , 

where Sij is the Kronecker symbol. 

3.1. Canonical differentials. We choose canonical holomorphic differentials (of the first kind) dit 4 = 
(dui, . . . , du g ) and associated meromorphic differentials (of the second kind) dr* = (dri, . . . , dr g ) in such 
a way that their g x g period matrices 



(3.4) 



a* 



2w = ( i> dui ) , 2oj' = ( &> du 



i,k=l,...,g V J b / i,k=l,...,g 



tlk 



2 V = ( - f dn) , 2<n' =(- & d n 

i,k=l,...,g V J[, fc /i,fe=l,...,g 



satisfy the generalized Legendre relation 



(3.5) MJA/ T = -yJ 
with 

with ff and l ff as the zero and unit g x (/-matrices. Such a basis of differentials can be realized as follows 
(see Baker (1897), p. 195) 

(3.7) du(z,w) = U ^ d \ U l (z)=z l -\ i = l...,g, 

w 

TZ(z)dz 2g+l-i 

(3.8) dr{z,w) = ^— , Ki(z)= ^ (fc + 1 - z)A fc+ i +l z fe , i = l..., ff , 

/c— i 

where the coefficients Xi are given by (|3.ip . 

We denote by Jac(A" s ) the Jacobian of the curve X g , i.e., the factor C ff /r, where r = 2lu ® 2w' is the 
lattice generated by the periods of the canonical holomorphic differentials. Any point u £ Jac(A g ) can be 
presented in the form 

(3.9) u = 2uje + 2Je' , 

where e, e' £ M. 9 . The vectors e and e' combine to a 2 x g matrix and form the characteristic of the point 
u, 



(3-10) [u] := , , r| 

If u is a half-period, then all entries of the characteristic e are equal to \ or 0. 
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Beside the canonical holomorphic differentials du we will also consider normalized liolomorphic differ- 
entials denned by 

(3.11) dv = (2uj)- 1 du. 

Their corresponding holomorphic periods are l g and r, where the Riemann period matrix r := u u 1 is 
in the Siegel upper half space & g of g x ^-matrices (or half space of degree g), 

(3.12) & g = |r g x g matrix|r T = r, Imrpositive definite} . 
The corresponding Jacobian is introduced as 

(3.13) lac(X g ) := (2^)^ Jac(A 3 ) = C 9 /l g © r . 

We will use both versions: the first one (2u>, 2a/) in the context of cr-functions, and the second one (l ff , t) 
in the case of ^-functions. 

The Abel map 21 : (X g ) n — > C 9 with the base point Po relates the set of points (Pi, . . . ,P n ) (which are 
called the divisor T>) with a point in the Jacobian Jac(X s ) 

" rPk 

(3.14) 2l(P,...,P„) := V / du. 

fe=l Jp o 

The divisor 2? in (|3.14[) can be also denoted as Pi + . . . + P n — rip). 
Analogously we define 

(3.15) 2l(P 1 ,...,P„)= V / d- y = (2 w )- 1 2l(Pi,...,P„). 

fe=l 7p o 

In the context of our consideration we take Po as infinity, Po = (oo,oo). 

3.2. 0— functions. The hyperelliptic ^-function with characteristics [e] is a mapping 8 : Jac(A g ) x & g — > C 
defined through the Fourier series 

(3.16) 9[e](v\r):= ^ ^{MM^l+aM*^)} _ 
It possesses the periodicity property 

(3.17) 9[e]{v + n + tu'\t) = e - 2i ™' f ( v +i Tn "> e 2i *( nte '- n ' te) 0[e](v\T) . 

For vanishing characteristic we abbreviate 0(v) := 9[0](v\t). 

In the following, the values e' k will either be or |. The equality p. 171) implies 

(3.18) 9{e}(-v\T)=e- 4 ™ t£ '8{e}(v\T), 

so that the function 0[e]('w|t) with characteristics [e] of only half-integers is even if 4e t e' is an even integer, 
and odd otherwise. Correspondingly, [e] is called even or odd, and among the 4 9 half-integer characteristics 
there are i(4 5 + 2 9 ) even and i(4 9 — 2 9 ) odd characteristics. 

The non- vanishing values of the ^-functions with half-integer characteristics and their derivatives are 
called (^-constants and are denoted as 



1[e] := 9[e](0;r), 9 h3 [e] ~ ^j^M^ T ) 

i[e] := ^-0[s}(z;r) 



, etc. for even [e]; 

z=0 



etc. for odd [e] . 



Z =0 ^*JlV*jL 

Even characteristics [e] are called nonsingular if 9[e] ^ 0, and odd characteristics [e] called nonsingular if 
9i[e] ^ at least for one index i. 

We identify each branch point e 3 of the curve X g with a vector 

(3.19) 2lj := / du =: 2w£j +2u'e' j G Jac(X 9 ), j = 1, . . . ,2g + 2 , 

J oo 

what defines the two vectors £j and e'-. Evidently, [2l2 9 +2] = [0] = 0. 

In terms of the 2^ + 2 characteristics [21 j] all 4 s half integer characteristics [e] can be constructed as 
follows. There is a one-to-one correspondence between these [e] and partitions of the set Q = {1, . . . , 2<7+2} 
of indices of the branch points ([Fay 73;, p. 13, [Bak897j p. 271). The partitions of interest are 

(3.20) l m U J m = {ii, . . . , i s+ i_2m} U {ji, . . . , j s +l+2m}, 
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where m is any integer between and [^-j^] ■ The corresponding characteristic [e m ] is defined by the 
vector 

S+l-2m 

(3.21) A m = J2 ^i k + K 00 =:e m + Te' m , 

k=l 

where Koo € 3ac(X g ) is the vector of Riemann constants with base point oo, which will always be used 
in the argument of the ^-functions, and which is given as a vector in Jac(X g ) by 

(3.22) K x := *i 

all odd [2lj] 

(see e.g. jFK80j . p. 305, for a proof). 

It can be seen that characteristics with even m are even, and with odd m are odd. There are | \g+i) 
different partitions with m = 0, (g 9 ^ ) different partitions with m = 1, and, in general, ( g _^f^2m) down to 
( 29 + 2 ) = 2g + 2 partitions if g is even and m = g/2, or ( 29 q 2 ) = 1 partitions if g is odd and m = (g+ l)/2. 
One may check that the total number of even (odd) characteristics is indeed 2 2s_1 ± 2 9_1 . According to 
the Riemann theorem on the zeros of ^-functions |Fay73| , 8(A m + v) vanishes to order m at v = Q and 
in particular, the function #(.R"oo + v) vanishes to order \ a ^\ at v = 0. 

Let us demonstrate, following FK80 , p. 303, how the set of characteristics [21^] = [2tfe], k = 1, . . . , 2g+2 
looks like in the homology basis shown in Figure 1. Using the notation f k = ^{8\k 1 ■ ■ ■ , Sgk) 1 and for 
the fc-th column vector of the matrix t, we find 



(3.23) 2l 29 =2l 2g+1 - | dv = ^f k 



1(0 ... 
2 1 1 ... 1 1 



e2 s + l fc ~ 1 



g 



r« i 1 /0 ... 1 
[^=2 1 1 ... 1 1 



3029-1 = 2t2 S - / - ^ / fe + r ff , -► [Oaa-ij 

Continuing in the same manner, we get for arbitrary 1 < fc < g 

k 



1/0 ... 1 
2 ll 1 ... 1 



[21: 



2fc+ 2 J 




1 

1 



(3.24) 



and finally 

o*> S ::: S), I ::: S 

The characteristics with even indices, corresponding to the branch points e 2 „, n = 1, . . . ,g, are odd (except 
for [2t 2g + 2 ] which is zero); the others are even. Therefore in the basis drawn in Figure [TJ we get 

g 

(3.26) ^ = ^2l 2fc . 

fc=i 

The formula (|3.26l) is in accordance with the classical theory where the vector of Riemann constants is 
defined as (see Fay |Fay73| , Eq. (14)) 

(3.27) Divisor K Po = A - (g - 1)P , 

where A is divisor of degree g — 1 that is the Riemann divisor. In the case considered Pq = oo and 
A = e 2 + e 4 + . . . + e 2g — oo. The calculation of the divisor of the differential FJ/Ui ( x ~ e 2k)dx/y leads to 
required conclusion 2A = JCx where K-x is canonical class. 
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3.3. a— functions. The Klcinian cr-function of the hyperelliptic curve X g is defined over the Jacobian 
Jac(X g ) as 

(3.28) cr(«; M) := C*6'[K oc ]((2cj)" 1 m; r)exp {u T xu) , 

where x = r/(2uj)~' 1 , the constant 

-1/4 

r= ,/ , " , I n ^ 

.l<i<j<2g+l 



(3.29) 



det(2w) 



and M defined in (I3.6[) contains the set of all moduli 2cj, 2uj' and 2rj : 2r/. In what follows we will use the 
shorter notation o~(u;M) = a(u). Sometimes the cr-function (|3.28|l is called fundamental cr-function. 

The multivariable cr-function (|3.28|) represents a natural generalization of the Weierstrafi cr-function 
given by 



(3.30) 



cr(u) 



2u ^(e 1 -e 2 )(e 1 -e 3 )( e 2-e 3 ) 1 6XP \ 2uj 



1 . 



where #j is the standard ^-function. 

The fundamental cr-function introduced by the formula p. 281) respects the following properties 

• it is an entire function on Jac(A" ff ), 

• it satisfies the two sets of functional equations 



(3.31) 



(3.32) 



a(u + 2uk + 2u/fc; M) = ^ivk+n'k'^u+uk+^k') ^ M ^ 

cr( U ;( 7 M T ) T ) = cr(M;M), 

where 7 € Sp(2p,Z), that is, jJj t — J, and M T is the matrix M with interchanged submatrices 
u)' and 77. The first of these equations displays the periodicity property, and the second one the 
modular property. 
• In the vicinity of the origin the power series of cr(u) is of the form 

cr(tt) = Stt(u) + higher order terms , 

where S n (u) are the Schur- Weierstrafi functions whose definition we will recall in the next sub- 
section. 



3.4. Schur— Weierstrafi functions. The Schur function is a polynomial in the variables u%,U2, ■ ■ ■ built 
by a partition A : Ai > A2 > . . . > A„ of weight |A| = X)"=i ^i- -^ n y partition A can be written in the 
Frobenius notation A = (011, . . . , a r \0i, . . . , /3 r ), where the number r is the rank of the partition and the 
integers (ctj,f3j) are the numbers of nodes in the Young diagram to the left from j-th diagonal node and 
down to it (see |Sag01| ). 

In what follows we will deal with the special kind of Schur functions that are related to the hyperel- 
liptic curve and that we will call Schur- Weierstrafi functions S n (u) following BEL99]. The associated 
partition ir is defined by the Weierstrafi gap sequence w = (wi, . . . ,w g ) at the infinite branch point, 
w = (1, 3, . . . , 2g — 1) by the formula 



(3.33) 



71* = Wg-k+i + k - g, i = L,...,g. 

In the considered case the associated Young diagrams are therefore symmetric and satisfy the constraint 
7Tfc — 7Tfc + i = 1, i.e., in the cases 5 = 5 and g — 6 we have the diagrams 



i = l,.. 



and 

corresponding to a partition of rank 3. In general, for arbitrary g the rank of the partition is the integer 
part[*±i]. 

For the polynomials S n (u) the following representation is valid 



(3.34) 



S n — det (c 7Ti —i+j)^ 
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where e k is given by the determinant 
(3.35) 



Pi 

P2 



1 
Pi 



Pk-1 Pk-2 Pk-3 
Pk Pk-l Pk-2 







jfe-1 
Pi 



and the quantities p k are related to the Jacobian variables (u%, . . . , u g ) as 



(3.36) 

For example, 



(3.37) 


9 


= 1 


(3.38) 


9 


= 2 


(3.39) 


f! 


= 3 


(3.40) 


f! 


= 4 



Pk = ku g _ [k/2 ], k=l,...,2g. 



Sl(ui) = Ml, 

S2,i(ui,u 2 ) = ^u 3 , - Ul, 



>S , 3,2,l( U l ; '"2,'«3) 
5 , 4,3,2,l(wl,U 2 ,M3,M4) 



45 



,10 



4725 



-U2 U 3 — u \ + u l U 3l 

m!m3 H — V? A Ui 

105 4 15 4 3 3 4 



+ u 2 u 3 ul — ul + UiU 3 . 



3.5. Kleinian p— functions. The p-functions are a natural generalization of the corresponding Weier- 
strafi functions and given as logarithmic derivatives of a 

d 2 

Pij(u) = - — ^ \na(u), 

(3.41) 



duiduj 

Q3 

pijk(u) = - a n a lncr(u) , etc. 



duidiijduk 

where k € {1, . . . , g}. In this notation the WeierstraB p-function is pi t %(u). For convenience, we will 
also denote the derivatives of the er-function by 



(3.42) 



d d 2 



dui 



duiduj 



etc. 



The Jacobi inversion problem is the problem of the inversion of the Abel map and can be formulated 
as follows: for an arbitrary vector u € Jac(X 9 ) find the symmetric functions of g points P\, . . . ,P g G X g 
from the equation (|3.14j) . that is u = Ylk=i f * du. In the considered case of a hyperelliptic curve Jacobi's 
inversion problem is written in coordinate notation as 



(3.43) 



r Pi dx 


dx 


JPc v 




f Pl xdx 


f Pa xdx 


Jp y 


J Pa y 


x 9_1 dx 

\-...-\ 


rP g x g-idx 


V 


Jpo y 



— = Ul 



u 2 ■ 



"g i 



where P k = (xkiUk)- It can be solved in terms of Kleinian p-functions, i.e., X\, 
solutions of 



, x g are given by the g 



(3.44) 



„<?-! 



and y k = -p ggg (u)x 9 k 



9-1 



p g ,i(u) =0, 
■~Pggl( u ), k = l,. 



Another expressions for symmetric functions Ylt—! x l k , with I = 1, 2, 3, . . ., are given in [Vanh95 . 

Among the various differential relations between the Kleinian p-functions we quote the representation 

_i_ 9 ( 9 ~f" 1 ) 

of the Jacobi variety as algebraic variety in C ffH 2 obtained in [BEL 97] and used in the foregoing 
development. It is described by the set of cubic relations that can also be represented as minors of a 
certain matrix |BEL97) . 



PggiPggk 



^PggPgiPgk ~ 2{pgiPg-l,h + PgkPg- 



^{pgkPg^-1 + PgiPg,k-l) 
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(3.45) +4pfc_ 1>i _ 1 - 2(pk,i-2 + Pi.k-2) + ^2gPgkPgi + ^| 1 {SigPgk + 0~kgPgi) 

+ A 2l -2<S 4 fc + ^(A2i-l<Jfc,i+l + A 2 fc-l<5j,fc+l), 1 < i, fc < 5 . 

For 3 = 1 Eq. (|3.45p reduces to the Weierstrafi cubic p' 2 — 4p 3 — g 2 p — 33 if we set A2 = 0, pm = 
p'j Pn = P- 

3.6. Period matrices. The construction of the fundamental a-functions requires the knowledge of the 
first period matrix (2w, 2a/) and the second period matrix (2r/, 2rj') which satisfy the generalized Lcgendre 
relation (|3.5[) . We will show that the second period matrix can be constructed in terms of the first period 
matrix and even ^-constants. 

For that purpose we consider any even nonsingular half-period of the ^-function (2w) 2lx + K x 
where the half-period 



(3.46) 2Ci = / du+...+ du 

J OG J OO 

corresponds to the partition of the branch points 

(3.47) {i 1 ,...,i g ,2g + 2}\J{j 1 ,...,j g+1 } =1 U J . 

Then from (|3.44p with Xk — e% kl u = 2lx , follows an expressions for pi ff (2lx ) in terms of symmetric 
functions of the elements , ik G Iq, 



(3.48) 



e»x + ■■■ + e lg = p 9S (2lx ), 
e^e^ + . . . + ei t&i = -p ff _i ifl (2li ), 



^•••ei, = {-i) 9 - l Plg {m Io ). 



One can see that all two-index symbols py(2lx ) can be expressed in terms of symmetric functions of 
two sets of variables e,, i G Iq and e^, j G J7b- That follows from the fundamental cubic relation (|3.45[) . 
Indeed pi flff (2lx ) = for all i = 1, . . . ,g and we get g ( g + 1 ) relations to which (|3.48[) are substituted as 
well as expressions for the X's in terms of e», i = 1, . . . , 2<7 + 1. These equations can always be solved with 
respect of the remaining pjfc (2lx ) . Therefore, all pz,j(2lxo) are known in terms of the branch points efe. 

Proposition 3.1. Let 2lx + 2i^K. OQ be an arbitrary even nonsingular half-period corresponding to the g 
branch points of the set of indices Iq — {ii, . . . , i g }. We define the symmetric 3x3 matrices 



(3-49) <p(2lxj :=(py(2lxo)) lJ=1 ,..., g 

and 

(3.50) T(2l Io ) := ("^^log *[*»](*; r)^-^ J 



TTien i/ie x-matrix is given by 

(3.51) x = -~9J(2lx ) - i((2c)- 1 ) T T(2lx )(2 W )- 1 

and i/ie half-periods r\ and r/' of the meromorphic differentials can be represented as 

(3.52) t]^2m:uj, rf = 2W - y (w^ 1 ) 7 . 

We remark that (|3.51|) represents the natural generalization of the Weierstrafi formulae 

2 1^'(0) _ n 2 i^'(o) „ _ „ 2 i<(0) 



(3.53) 2^ = -2e^--^, 2 w = -2e 2 ^--^, 2^ = -2e^ 2 ^ , 

see e.g. the Weierstrafi-Schwarz lectures, Wci893 p. 44. Therefore Proposition |3~T1 allows the reduction 
of the variety of moduli necessary for the calculation of the a- and p-functions to the first period matrix. 
The generalization of the result fl3.5ip to non-hyperelliptic curves is considered in [EKllj . 

4. Inversion of one hyperelliptic integral 



Classically it is known that only symmetric functions of 3 points of algebraic curves of genus 3 can be 
presented as single-valued functions over the Jacobi variety of the curve. Nevertheless, one hyperelliptic 
integral can also be inverted analytically in terms of c-functions restricted to the 0-divisor. 
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g 


m(0 o ) 


m(6i) 


m(§ 2 ) 


m(0 3 ) 


m(§ 4 ) 


m(0 5 ) 


m(9 6 ) 


1 


1 















2 


1 


1 





- 


- 


- 


- 


3 


2 


1 


1 











4 


2 


2 


1 


1 









5 


3 


2 


2 


1 


1 







6 


3 


3 


2 


2 


1 


1 






Table 1. Orders m(0fc) of zeros 9(Qk + v) at v = on the strata 0fc. 



4.1. Stratification of the 9— divisor. The ^-divisor 9 is defined as the subset of J&c(X g ) that nullifies 
9 and, therefore, the er-function, i.e. 



(4.1) 



e 



{. 



€ Jac(X g 



The subset 0^ C O, < fc < g, is called fc-th stratum if each point »£6 admits a parametrization 



(4.2) 



dt> + if 



where ©o = {-RToo} and g -i = 0. We furthermore denote 9 = Jac(X g ). The following natural 
embedding is valid: 

(4.3) © C 6i C . . . C ff _i C Gg = 3^C(Xg) . 

We define the ^-function to be vanishing to the order m(Ok) along the stratum 0fc if for all sets otj, 
j = 1, . . . , g with 0<ai + ... + a g < m 

gai + ...+a g 



(4.4) 



. di 



while there is a certain set of a, , with <x\ 



(v\t) = 0, V«€6 t 



= mso that (|4.4[) does not hold. The orders m(0j,) 



of the vanishing of 9(0^ + v) along the stratum 0^ for the first genera are given in Table [TJ 

In the following we focus on the stratum Qi corresponding to the variety Qi C Ja.c(X g ), which is the 
image of the curve inside the Jacobian, 



(4.5) 



©i 





f p } 


L G 


u = du 




J OG J 



and 0i = (2a;)" 1 ©! +K C 



We remark that another stratification was introduced in Vanh95 for hyperelliptic curves of even order 
with two infinite points oo + and oo_ that was implemented for studying the poles of function on Jacobians 
of these curves. The same problem relevant to strata of the 0-divisor was studied in [AFOOj . 

4.2. Inversion formulae. For the case of genus two the inversion of a holomorphic hyperelliptic integral 
by the method of restriction to the 0-divisor was obtained independently by Grant [Gra90| and Jorgenson 
|Jor92j in the form 



(4.6) 



o-i (u) 



(72 (u) 



u = (u 1 ,u 2 ) 1 



<t(u)=0 



This result was implemented in [EPR03 , and explicitely worked out in the series of publications HL08 , 
IHL08a[ IHKKL081 IHKKL091 lHKKL09aj . and oth ers. 

The case of genus three was studied by Onishi Oni98 , where the inversion formula is given in the form 

o-ia (u) 



(4.7) 



023 (") 



, U = (ui,U2 1 U 3 ) T . 

a(u)=<j 3 (u)=0 

Formula (|4.7I) is based on the detailed analysis of the genus three KdV hierarchy and its restriction to the 
6*-divisor. Below we will present the generalization of (|4.6|) and (|4.7|) to higher genera. For doing so we 
first analyze the Schur-Weierstrafi polynomials that represent the first term of the expansion of o~(u) in 
the vicinity of the origin u ~ 0. 
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4.3. Schur Weierstrafi polynomials on 6^-strata. The Schur-Weierstrafi polynomials associated to 
the curve X g are defined in the space <C 9 9 (ui, . . . , u g ) by the Weierstrafi gap sequence at the infinite 
branch point that for g > 1 is always a Weierstrafi point. The 0-divisor O and its strata ®k in the vicinity 
of the origin u ~ are given as polynomials in u. An analysis of the Schur-Weierstrafi polynomials leads 
to the following 

Proposition 4.1. The following statements are valid for the Schur-Weierstrafi polynomials S n (u) asso- 
ciated with a partition 7T G 

(1) In the vicinity of the origin, an element u of the first stratum 8i C 6 is singled out by 

d j 

(4.8) S„(u) = 0, _ T 5 W («) = :j 1 g 2. 

OUg 

(2) The derivatives fulfill 



d J f = if 1 < i < SiS^H 

(4.9) ^ s "<"»{; l " EGl 

(3) The following equalities are valid for u G 0i 

QM qM+1 
{ ' ul~ d M , , ~ d M+1 



where M = \{g - 2) (5 - 3) + 1. 
(4) The order of vanishing of restricted to &% is the rank of the partition tt. 

It was noted in [BEL 9 9] that the Schur-Weierstrafi polynomials respect all statements of the Riemann 
singularity theorem. In particular, if 



(4.11) Z=\ ...,—,z 

V ' \2g- 1' '2fc-l ' 3 ' 

and if 7r is the partition at the infinite Weierstrafi point of the hyperelliptic curve X g of genus g, then the 
function 

(4.12) G(z) := S n {Z-u) 

either has g zeros or vanishes identically. This result was extended in jMP08 ( MP10J. Moreover we will 
conjecture here that the properties of the Schur-Weierstrafi polynomials given in Proposition 14.11 can be 
"lifted" to the fundamental tr-function (|3.28p . 



4.4. Inversion for higher genera. The above analysis permits to conjecture the following inversion 
formula for the general case of hyperelliptic curves of genus g > 2 



d 



M+l 



(4.13) 

and 
(4.14) 



dmdu^ 



cr{u) 



qm+i 
du-iduf 



M 



(9-2)(g-3) 



+ 1 



ue&i 



81= tie Jac(X g ) a(u) = 0, -a(u) = V j = 1, 

I du}, 



The analog of this formula for strata 0^, 1 < k < g and (n, s)-curves in the terminology of (BEL99 
was recently considered by Matsutani and Previato [MP08j . [MPlOj . 

We remark that the half-periods associated with one of the branch points e\, . . . , &2g+\ are elements of 
the first stratum and, therefore, the following proposition is valid: 
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Proposition 4.2. Let the curve X g be of genus g > 2 and let 21^ be the half-period that is the Abelian 
image with the base point Pq = (oo,oo) of a branch point ej. Then 

qm+i 

(4.15) e t = , 

where M = \{g - 2)(g - 3) + 1. 

The formula (|4.15|) can be considered as an equivalent of the Thomae formulae |Tho8 70 . Similar 
formulae can be written on other strata 0^, but we only mention here the set of formulae that we are 
using in our construction. 

Proposition 4.3. Let X g be a hyperelliptic curve of genus g and consider a partition 

(4.16) IrU Ji ={i 1) ...,i 9 _i}U{ii,...,i 9+3 } 
of branch points such that the half-periods 

(4.17) (2w)- 1 2l Il +KOO e e H U8 r2 
are nonsingular odd half-periods. Consider two cases: 

(4.18) I[ :={i 1 ,...,i g - 1 }J2g + 2, 

(4.19) X^':={^ 1 ,..., ^^ _ 1 }^2 5 + 2, i.e. t g - 1 =2g + 2. 

Denote by Sfe(X 1 ) and Sfc(X") the elementary symmetric function of order k built by the branch points 
ejj , . . . , &i ^ and , . . . , ei g _ 2 correspondingly. Then the following formulae are valid 

s k (l[) = (-i)fc^l , k = l,...,g-l, 

(4-20) a l . 

Sfc _ 1 (^) = (-l)^^zi(2l z ,) , k = 2,...,g-l, 

where we denoted for typographic convenience ^ ^ ='■ ^(2l)- 
The following corollary follows immediately from (|4.20|) : 

Corollary 4.4. Let g > 3 and I' = I" U {i} where I" — {ii, . . . , ig-2} and i ^ 2g + 2, i ^ I". TTien £/ie 
representation for the branch points ej 

(4.21) ei = _^zi( 2ll;) + ^zi( 2ll „). 
is valid. 

The formulae f|4.21[) can be understood as a generalization of those given in Bolza B0I886. We also 
remark that a comparison of the two representations (|4.15[l and (14.21[) of the branch point leads to an 
interesting 0-constant relation. 

4.5. Calculation of moduli. Calculations in terms of 9- or er-functions are usually considered as tech- 
nically cumbersome what prevents wide applications of algebro-geometric methods. In particular, the 
procedure of the evaluation of the period matrix in the given homology basis is technically complicated 
even in the case of a hyperelliptic curve. Based on the above analysis we show that modern software like 
the "Maple/ algcurves" package now allows the calculation of the 9- or a-functions without drawing and 
even without knowledge of the homology basis, at least in the hyperelliptic case. The calculation scheme 
is given by the following steps: 

Step 1.: For the given curve compute first the period matrices (2u>, 2u/) and r = w _1 u/ by means of 
the "Maple/algcurves" code. Compute then the winding vectors, i.e., the columns of the inverse 
matrix 

(4.22) {2U)- 1 = {U ll ...,U g ). 

Step 2.: We then find all nonsingular odd characteristics by direct computation of all odd 9— 
constants. According to Table Q] we have two sets B\ C 9 -i and B2 C 9 -2 of nonsingular 
odd half-periods. For each element of b\ € B\ there are , . . . , e% t =/= 00 such that 

(4.23) 61 = [ n dv + ...+ f '"^ dv + K^e Q g -! 
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and for each element of 62 € B% there are e» 



7^ 00 such that 



(4.24) 



dv 



dv + G 9„_2 



By using (|4.20p and the known values of the winding vectors one can find the correspondence 



between the sets {e^,...,^ t } and {f 



,} of branch points and the nonsingular odd 



characteristics [(2a 



2L 



Koo] and [(2a;)- 1 (2l 4 



Then one can add 



(4.25) 



these characteristics and find the one-to-one correspondence 

" ! d^ [2L J, i=l,...,2$- 



Step 3.: Among the 2g + 2 characteristics (|4.25l) there should be precisely g odd and g + 2 even 
characteristics. The sum of all odd characteristics gives the vector of Riemann constants with the 
base point at infinity. Check that this characteristic is of order ■ 

Step 4.: Calculate the symmetric matrix x by (|3.5ip and then the second period matrices 2rj, 2r/ 
according to the Proposition 13. II 

We add two remarks: 

Remark 1.: The proposed way to compute the [21&] is not the only possible way. One can also 
use the standard Thomae formulae for the ^-constants with even characteristics while we used 
Bolza-type formulae for odd characteristics. Also the "Maple/ algcurves" contains expressions for 
homology cycles written in terms of paths connecting branch points and it could be possible in 
principle to find the [21&] solving the system of equations. 

Remark 2.: ] One could imagine a case where the drawing of the homology basis is nevertheless 
important for the problem, e.g., in the case when the curve possesses additional symmetry. In this 
case one can use T. Northower's program [Nor 101 INorlOaj to find the symplectic transformation 
between the required basis and the one given by "Maple/ algcurves" . 

4.6. Procedure of the inversion. We are now in the position to calculate the inversion, that is x(t), of 
the hyperelliptic integral 



(4.26) 



z k dz 



t. 



0<k<g. 



This will be carried through along the following steps. 

Step 1.: We first fix the homology basis, e.g., the basis of "Maple/algcurves" and compute all moduli 

of the curve according to Sec. 14.51 
Step 2.: The dynamic system considered is evaluated between two branch points of the polynomial 

*P2g+\{%) that defines the curve (|3.1[) . Therefore fix a branch point, say ej, that is the starting 

point of the system evolution and find the half-period 21^ = J e * du. 
Step 3.: Use the formula (|4.15j) for g > 2 or (|4.6[) for g = 2 with the argument u of the cr-function 

given as 



(4.27) 



u = 2L 



/fc-iW 
t 

fk+l(t) 



with /(0) = . 



V /,(*) / 

where f(t) = (fi(t), . . . ,f g (t)) T are locally given functions that resolve the conditions (|4.14l) of 
the restriction to the stratum ©x- The vector function f(t) can be obtained from /(0) using the 
Newton method. In this case the approximation process should be carried through for the real 
and imaginary parts of each function fi(t) separately. 

We emphasize that the inversion procedure given above carries local character and the quasi-elliptic 
function can be defined only locally. But we believe that our method elucidates the geometric structure of 
the object and leads to exact calculations in contrast to numerically solving ordinary differential equations. 

Below we will consider the case g = 2 and g = 3 and demonstrate in more detail how this scheme can 
be applied. 
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Figure 2. Homology basis on the Riemann surface of the curve X 2 with real branch 
points e\ < e 2 < . . . < e$ = 00 (upper sheet). The cuts are drawn from e2,:_i to e 2 i, 
i = 1,2,3. The b-cycles are completed on the lower sheet (dotted lines). 



5. Hyperelliptic curve of genus two 
We consider a hyperelliptic curve X 2 of genus two 

w 2 = 4(z - ei)(z - e 2 )(z - e 3 )(z - e 4 )(z - e 5 ) 
= 4z 5 + A 4 z 4 + A 3 z 3 + X 2 z 2 + X 1 z + A . 
From p.7p and (|3.8j) the basic holomorphic and meromorphic differentials are 



(5.1) 



(5.2) 
(5.3) 



dz 

dui = — , 

w 

zdz 

du 2 = , 

w 



12z 3 + 2A 4 z 2 + X 3 z , 
dri = : d^ . 



4u; 



Then the Jacobi inversion problem for the equations 



(5.4) 

is solved in the form 
(5.5) 



iz 

w 

{zi ' Wl) zdz 
w 



dr 2 = — dz . 

w 



(z 2 ,w 2 ) dz 

— = ui , 

w 



00 



zdz 
w 



u 2 



z\ + z 2 = p 22 (u) 7 ZiZ 2 = -pi 2 (w) , 

w k = -p222(u)z k - p 122 (u), k = l,2. 



5.1. Characteristics in genus two. The homology basis of the curve is fixed by defining the set of 
half-periods corresponding to the branch points. The characteristics of the Abelian images of the branch 
points are defined as 



(5.6) 

what can be also written as 



du 



2l t = 2uj £i + Zoj'e'i, 
In the homology basis given in Figure [5] we have 



f' f' 



1,...,6. 



(5.7) 
(5.8) 



[SZti] = 



1 (\ 



2 10 



[2l 2 ] = 



1 A 



[sa 4 l = :r 



1 /o 1 



2 V 10 
1 /o 



2 VI V [ * 5]= 2{l 1 

The characteristics of the vector of Riemann constants Koo is 

(5-9) [«oo] = [St 3 ] + [SU] = 5(j J 



[sfe] = 

[«6] = 



1 /o 1 

2\l 

1 Ad 
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From the above characteristics 16 half-periods can be build as follows: Denote the 10 half-periods for 
i ^ j = 1, . . . , 5 that are images of two branch points as 

(5.10) Sl i , j = 2u(e i + e j ) + 2u)'(e' i + e' j ), i = l,...,5. 
Then the characteristics of the 6 half-periods 

(5.11) [(2w)- 1 2t 1 + K x ] =:S U i = 1, . . . , 6 
are nonsingular and odd, whereas the characteristics of the 10 half-periods 

(5.12) [(2w)- 1 a j + K^} e itj , 1 < i < j < 5 
are nonsingular and even. 

Odd characteristics correspond to partitions {6} U {1, . . . , 5} and {k} U . . . , «4, 6} for i%, . . . , in ^ k. 
The first partition from these two corresponds to 0o and the second to 0i. 

From the solution of the Jacobi inversion problem we obtain for any i, j = 1, . . . , 5, i ^ j 

(5.13) e { + ej = p22(^ij), -eiej = pi 2 (l\j) . 
Using the relation (see }Bak903j . |BEL97j ) 

(5-14) p222 = 4 P22 + 4pl2p22 + 4pn + A3P22 + A4P22 + A 2 

one can also find 

(5.15) eie 3 (e p + e q + e r ) + e p e q e r = pn(ilij) , 

where i,j,p,q, and r are mutually different. 

From (f57T3|) and (jSTSl we obtain an expression for the matrix h that is useful for numeric calculations 
because it reduces the second period matrix to an expression in the first period matrix and ^-derivatives, 
namely, in the case et = e±, Bj = e2, 

/,,m l ( e 1 e 2 (e 3 + e 4 + e 5 ) + e 3 e 4 e 5 -e x e 2 \ 1 f0 wi T o-/r> Vo \-l 

(5-16) x=--\ ' ^ o )--(2u}) 1(n lt 2){2w) , 



-eie 2 ei + e 2 

where 1 is the 2 x 2-matrix defined in Proposition 13.11 

5.2. Inversion of a holomorphic integral. Taking the limit z 2 — > 00 in the Jacobi inversion problem 
(|5.4[) we obtain 

f {z ' w) dz f {z ' w) zdz 
5.17 / — =«i, / =u 2 . 

Jog W Jog W 

The same limit in the ratio 

pi 2 (w) z x z 2 



(5.18) 

P22(U) Zi+Z 2 

leads to the Grant- Jorgenson formula (|4.6[) . In terms of ^-functions this can be given the form 
(5.19) z dueiK^m^u-r) 



dvOlK^^-^T) 



J((2w)- 1 m:t)=0 



where here and below djj = X)f=i ^i~§r 1S ^ ne derivative along the direction U. Here we introduced the 
"winding vectors" U, V as column vectors of the inverse matrix 

(5.20) (2w) _1 = (U, V) . 

From (|4.6[) we obtain for all finite branch points 



(5.21) * = -^75n. < = 1 5 

or, equivalently, 

f5 22) _du0\M 

(5.22) e t - t -l,...,5. 

This formula was mentioned by Bolza [B0I886 (see his Eq. (6)) for the case of genus two curve with finite 
branch points. 
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e 8 = oo 



Figure 3. Homology basis on the Riemann surface of the curve X 3 with real branch 
points ei < e2 < . . . < e 8 = 00 (upper sheet). The cuts are drawn from e 2 i-\ to e 2i , 
i = 1,2,4. The b-cycles are completed on the lower sheet (dotted lines). 



6. Hyperelliptic curve of genus three 

As the next case we consider the hyperelliptic curve X 3 of genus three with seven real zeros as a model 
problem in advance to the real physical problems studied in the next section. Let the curve X 3 be given 

by 



(6.1) 



w 2 = 4(z - ei)(z - e 2 )(z - e 3 )(z - e 4 )(z - e 5 )(z - e 6 )(z - e 7 ) 



= 4z 7 + A 6 z 6 + ... + Aiz + A . 
The complete set of holomorphic and meromorphic differentials with a unique pole at infinity is 



(6.2) 



dz 

du! = — , 

w 

zdz 

du-2 = , 

w 

z 2 dz 
du 3 = 



dz 

dn = z(20z 4 + 4A 6 z 3 + 3A 5 z 2 + 2A 4 z + A 3 ) — , 



4w 



dz 

dr 2 = z 2 {Ylz 1 + 2X 6 z + A 5 ) — 

4w 



dr 3 = 



z 3 dz 



w 



Again we introduce the winding vectors 

(6.3) (2LO)- 1 = (U,V,W) . 

The Jacobi inversion problem for the equations 

rZl dz f Z2 dz f 23 dz 



w 



(6.4) 



is solved by 



f Zl zdz f Z2 zdz f Z3 zdz 

/ + / + / =u 2 , 

Joo W W Joe W 

r zi z^dz r Z2 z^dz r Z3 2 ' 

J 00 W Joo w J 00 



= u 3 



(6.5) 



zx + z 2 + z 3 = p 33 (u), z 1 z 2 + z 1 z 3 + z 2 z 3 = -p 23 (u), ziz 2 z 3 = p 13 (u) 
w k = -p 333 (u)zl - p 233 (u)z k - pi 33 {u), k = 1, 2, 3 . 



6.1. Characteristics in genus three. Let be the Abelian image of the fc-th branch point, namely 
(6.6) 2t fe = / du = 2uje k + 2uj'e' k , fc = l,...,8, 

J OO 

where e& and e f k are column vectors whose entries Skj, s f k j are \ or zero for all k = 1, . . . , 8, j = 1, 2, 3. 
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The correspondence between branch points and characteristics in the fixed homology basis is given as 

s xv w-if; » »v «=v o 1 



The vector of Riemann constants .K^ with the base point at infinity is given in the above basis by the 
even singular characteristics, 

(6-8) [Koa] = [^] + m + m = l(l 1 

From the above characteristics 64 half-periods can be built as follows. Start with singular even charac- 
teristics, there should be only one such characteristic that corresponds to the vector of Riemann constants 
Koo. The corresponding partition reads I2UJ2 = {}U{1, 2, . . . , 8} and the ^-function ^(.Koo +v) vanishes 
at the origin v = to the order m = 2. 

The half-periods Ai = (2w) _1 2lfc + K^ G Oi correspond to partitions 

(6.9) IxUJt ={fc,8}U{ii,...,j 4 }, ji,... J* i {8,fc} 

and the ^-function 0(Ai + v) vanishes at the origin v = to the order m = 1. 
Also denote the 21 half-periods that are images of two branch points 

(6.10) n i>j = 2u(e i + e j ) + 2u'(e' i + e' j ), i,j 1 7./../. 

The half-periods Ai = (2ui)~ 1 flij + -Koo G 02 correspond to the partitions 
(6.H) XiU Ji = {i,j}U{ji,...,H}, h,.-.,34<t{i,j} 

and the ^-function 0(Ai + v) vanishes at the origin, v = 0, as before to the order m = 1. Therefore the 
characteristics of the 7 half-periods 

(6.12) [(2w)- 1 2l 4 + fs: oo ] =:<$;, i = l,...,7 
are nonsingular and odd as well as the characteristics of the 21 half-periods 

(6.13) [{2w)- l n id + Kooj =: Sij , 1 < i < j < 7. 
We finally introduce the 35 half-periods that are images of three branch points 

(6.14) il i>jjk = 2uj(£i + ej + e k ) + 2uj'(e'i + e'j + e' k ) G Jac(A g ), 1 < i < j < k < 7 . 

The half-periods Ai = (2ui)~ 1 flij t k + Koo correspond to the partitions 

(6.15) l Uj = {i,j,k,8}U{j 1 ,...,j 4 }, j u . . . ,j 4 £ {i,j, k, 8} . 

The ^-function 6*(Ai + v) does not vanish at the origin v = 0. 
Furthermore, the 35 characteristics 

(6.16) [e i<jtk ] = [(2w)- 1 n iJ)fc + K^] , 1 < i < j < k < 7 

are even and nonsingular while the characteristic [.Koo] is even and singular. Altogether we got all 64 = 4 3 
characteristics classified by the partitions of the branch points. 

6.2. Inversion of a holomorphic integral. All three holomorphic integrals, 

f*dz [ x zdz [ x z 2 dz 

(6.17) / — =ui, / = tt 2 , / =u 3 



are inverted by the same formula (|4.7|) . Nevertheless, there are three different cases for which one of the 
variables Hi, 1*2,1*3 is considered as independent while the remaining two result from solving the divisor 
conditions <r(u) = (73(11) = 0. 

Formula (|4.7[) can be rewritten in terms of ^-functions as 

% ^[K^m^u) + 2(d u e{K 00 ]((2uj)-^u))e^u 
[ ' ' ' d^ w 9[K oc }((2u)-^u)+2(d v e[K <x }((2uj)-^u))e^Ku' 

where e 3 = (0, 0, 1) T . This represents the solution of the inversion problem. 
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From the solution of the Jacobi inversion problem follows for any l<i<j<k<7, 

(6.19) ei + ej + e k = p33(^ij,fc), —e.%ej — e^k - e 3 -ek = p23(^i,j,k), e.i^k = Pi3(&i,j,k) ■ 
From (|3.45j) we get the relations 

P333 = 4 P33 + ^6P33 + 4 P23P33 + A 5 p 33 + 4p 2 2 ~ 4p i3 + A 4 , 

(6.20) p| 33 = 4p^ 3 p 33 + A 6 p2 3 - 4p 22 p23 + 8pi 3 p23 + 4pn + A 2 , 
Pl33 = 4 P13P33 + A 6 p 13 - 4pi2pi3 + A 

and also find 

(6.21) Pi2(^i,j,fc) = S3S1 - S4, pu(£li,j,k) = S3S2 + S1S4 , p22(fii,j,fe) = S 3 + 2s 3 + s 2 Si , 

where si are the elementary symmetric functions of order I of the branch points ei,ej,ek and Si are the 
elementary symmetric functions of order I of the remaining branch points {1, . . . , 7} \ {i, j, fc}. 
From (|6. 19|) and (|6.21|) one can find the following expression for the matrix h 

(6-22) x = -iqj((fi iJ>fe )) - l -(2Lo)- lT %{^, h k ){2u)- 1 

where i,j, k are arbitrary and ^(fljj^) is the 3 x 3-matrix defined in Proposition 13. 11 

(6.23) i(n w ,n)= (- 92 i og 0[js: oo ]((2^)- 1 a J ,fc;T) 



m.n— 1,2,3 



For the branch points ei, . . . , eg the expression 
(6.24) e, 



du 


d w 


4- 22lf xe 3 




d v 


d w 


+ 22lf xe 3 


0[K oo ]((2 W )- 1 2l i ;r) 


J = 


1,.. 






&i + 


e 3 = 


^(^ij) 
0-3(0^) 




e, 


e i = 


O'l(^ij) 


. du8[S hJ ] 
' dw8[8ij] ' 






(7l(2li) 


du,0[Si] 


e 


i ~ 


<7 2 (Oi) 


d v 0[Si] ' 



(6.25) 

and for i = 1, . . . , 7 
(6.26) 

Using these data one can reconstruct by the known matrices 2w,2o/ all the characteristics [9tj], j = 
1, . . . , 2<? + 2, as well as the vector of Riemann constants following the procedure given in the Sec. 14.51 For 
example, for the genus 3 curve 

(6.27) y 2 = Ax{x - l){x -2){x- 3) (a: - 4){x - 5)(i - 6) 

we construct the "Maple/ alcurves" homology basis and compute the Riemann period matrices (2w, 2u/). 
We order the branch points according to 

(6.28) e\ = 0, e-2 = 1, e 3 = 2, e 4 = 3, e 5 = 4, e 6 = 5, e 7 = 6, e 8 = 00 . 
Then we find the Abelian images of the branch points 



2 V 1 1 v ' 2 V 1 V ' 2V 001 

There are 3 odd among these characteristics and therefore 

1/111 



(6-30) [K^ } V| (| | 

One can check that [K^] is the only even characteristic such that O^K^,]) = as follows from Table [T] 
Now one can construct arbitrary nonsingular even characteristics and calculate x and then the second 
period matrix. 
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7. An application: 9-dimensional Reissner-Nordstrom-de Sitter space-time 

Being armed with the method developed we now come back to one of the physical model systems 
presented in Section [2j We consider here the geodesic equation in a Reissner-Norstrom-de Sitter space- 
time of nine dimension. Setting d — 9 and introducing a new dimensionless coordinate f = r/rg and the 
dimensionless parameters L — L/rg, q — q/rs, and A = Arg we obtain from (|2.11j) 



(7.1) 
where 



df_y _ R 16 (r) 
dip / 



(7.2) L 2 R l6 {f) = f 16 S^- + \E 2 -S + ^-L 2 ) r 1A - r 12 L 2 + r 8 6 + f 6 L 2 - q 2 5f 2 - q 2 L 2 . 

28 \ 28 I 

Through the substitution f — l/y/u we halve the order of the polynomial and by u — x _1 + us, where us 
is any root of the polynomial Rg(u), we obtain a polynomial of seventh order so that Eq. (|7.ip reduces to 

/ d \ 2 7 7 

(7.3) \ x d) ='Pr(x) = ^2b i x i = brl[(x-e i ). 

^ ^ ' i=0 i=l 

Thus, solving the differential equation is reduced to the inversion of a genus three holomorphic hyperelliptic 
integral 

(7.4) Ifi - ifi in = 



where ipi n denotes the initial conditions for if and x- m is the starting point of the integration. 
From (|4.7j) we find the solution of the equation of motion (|7.1j) 



(7.5) r(f) 



Q~23(m) 
0"13(«) 



- u 8 




where 

(7.6) w = 2l s: +l /a(^-Wn) I, /i(0) = / 2 (0) = 

\ <P - Vin J 

and where the functions fi(f — fin) and ]i(ip — <f\-a) can be found from the conditions a(u) = and 
173(14) = 0. Also Xin is chosen as a branch point of the polynomial Vt{x) which defines the half-integer 
characteristic 21^ . The homology basis is fixed by the characteristics (|6.7[) and the cr-quotient is computed 
according to the formula (|6.18[) while the genus three ^-functions were calculated with the aid of the 
"Maple/RiemannTheta" code. 

The structure of the orbits is given by the number of zeros of the polynomial V7. This depends on 
the choice of the parameters E, L, A, and q. We choose a certain nontrivial value for A and q and 
draw the E 2 — L 2 parameter plot presented in Fig. | ^a)[ In the dark area the polynomial V7 possesses 
3 positive zeros resulting in one bound and one escape orbit, and in the brighter area it possesses one 
positive zero, yielding one two-world escape orbit. For certain values for E 2 and L 2 we obtain analytically 
the orbits shown in Fig. |4] From the causal structure of the metric encoded in the zeros of gtt = 5V?. 1 one 
concludes that after traversing both horizons a second time the test body enters a new universe (see, e.g., 
the discussion in jHKKLQ8j ). 

Furthermore, the observable perihelion shift ri pe riheiion can be given by a complete hyperelliptic integral 

(7.7) ^perihelion = 2 / / , A V ~ 2lT ■ 

Jr min y/V 7 (y) 



The perihelion shift compares the period of the radial motion with 27T. 
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(a) _E 2 -Z 2 -diagram: the dark gray re- 
gion indicates many— world peri- 
odic bound and escape orbits, as 
e.g. the orbits [(b)] and [(cj| In the 
light gray region there is one two- 
world escape orbit. 




(b) many— world periodic bound orbit 
(for a review of the classification 
of the geodesies in Reissner— 
Nordstrom— (anti) de Sitter 
space— times see 1HKKL08] .) 




(c) escape orbit 



Figure 4. Examples of test particle motion in a 9-dimensional Reissner-Nordstrom-de 
Sitter space-time. The chosen parameters for the cosmological constant and the charge 
are A = 8.7T0~ 5 , q = 0.4. We also consider massive particles, 5 = 1. For these parameters 
(a) shows the E 2 — L 2 parameter plot where the dark area indicates three positive real 
zeros and the bright area 1 real positive zero of Vj. For the chosen values E 2 = 1.045 and 
L = 0.5 two orbits are calculated analytically and plotted in|(b)|and (c) The black circles 



are the event and the Cauchy horizons. In this example e\ < &2 < e-a < e-A < Re(es) < ej, 
and es = ~e%. For the many- world periodic bound orbit [(b)] we choose a starting point e%, 
which corresponds to the half period [53jj] in (|6.7p . and the motion is bounded by e\ and 
e2- For the escape orbit (c) the starting point is which corresponds to the half period 
[2I3] in (|6.7p . and the test particle moves to infinity. 



8. Conclusion 

In this paper we presented the analytical solution of geodesic equations in general relativistic models, 
containing a hyperelliptic curve of arbitrary genus. Based on the solution for genus 2 and genus 3 and 
the investigation of the Schur-Weierstrafi polynomials, we present the solution (|4.15p for the geodesic 
equations with an underlying curve of arbitrary genus which is of the form or can be reduced to the form 

f x 

(8-1) ip-(p in = y l =, i = 0,...,g-l, 



y/P2g+i(y) 

where g is the genus of the curve w 2 = V2 g +i(x). As an example we integrate the geodesic equations 
with an underlying hyperelliptic curve of genus 3 in the Reissner-Nordstrom-de Sitter space-time in 9 
dimensions. 

One major task in the calculation of the analytical solution is the calculation of the first and second 
period matrices. In the Proposition 13.11 we provide a convenient and quick method for the calculation 
of the matrix x and the second period matrix from the first period matrix. The matrix h appears in 
the cr-function and, thus, in the general solution (|4.15[) . In this method no integration of meromorphic 
differentials is required and the result is given in terms of the holomorphic period matrix and f?-constants. 

In Section 14.51 we also propose a step by step algorithm for the calculation of the characteristics of the 
branch points and the vector of Riemann constants which appear in the 0— and c-functions and which 
are required for the calculation of the matrix x. We note that all the calculations can be done with 
the "Maple/algcurves" package, working in the homology basis automatically chosen by the computer 
program. 

The proposed solution (|4.15j) for the curves of g > 3 together with the solution for genus 2 curves has a 
wide spectrum of applications. These range from the geodesic equations in a wide class of black holes space- 



times like Schwarzschild-de Sitter |HL081 lHL08a) . Kerr-de Sitter [HKKL09aj . NUT-de Sitter |NdSprep 
and general type D Plebanski-Demianski HKKL09 space-times in 4 dimensions to the higher dimen- 
sional spherically symmetric Schwarzschild-(anti-)de Sitter, Reissner-Nordstr6m-(anti-)de Sitter space- 
times HKKL08 , higher dimensional axially symmetric Myers-Perry space-times |MPprep] discussed in 
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Section [21 and general Kerr-NUT-(anti-)de Sitter metrics in all dimensions introduced in [ChLP06 . This 
might be extended to space-times with cosmic strings HHLS10 , even in higher dimensions. The mathe- 
matical methods described here are also applicable to such problems as the motion of test particles with 
spin jHLSlOj and mass multipoles and to the motion of test particles with and without spin and mass 
multipoles in gravitating mass multipole fields, which have applications in astrophysics, satellite dynam- 
ics, and geodesy. The motion of test particles in the gravitational field of a black hole with disturbances 
related to mass multipoles have been discussed as test of the no-hair theorem (W09 . As a consequence, 
there is a wide range of applications of the formalism developed in the article to problems in the area of 
General Relativity. 
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